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PLASMAPAUSE AND DIFFUSION 


ACROSS FIELD LINES 


INTRODUCTION 

Whistler data (Carpenter 1963, 1966) and ion spectrometer measurements 
(Taylor et al., 1965) have shown, that the electron density decreases abruptly and 
very rapidly above L values of 3 to 5. Recent observations from OGO (Taylor, 
private communication) confirm this feature of the magnetosphere which is re- 
ferred to as the plasmapause. 

In a paper by Mayr, Brace and Dunham (1967), a study was conducted with the 
aim of describing and understanding the latitudinal variation of the ion composi- 
tion (Thomas et al., 1966, Taylor et al., private communication) for the two major 
ion constituents 0 + and H + during solar minimum. It was found that up to 45° 
geomagnetic latitude, the latitudinal increase of the electron temperature T e 
(Brace, Reddy, Mayr 1967) can properly account for the increase of the ion com- 
position transition level. Towards higher latitudes, however, the observed ion 
composition shows a rapid increase of the transition level which was in contrast 
to the theoretical result, which indicated a slight decrease of this level. In order 
to understand this discrepancy Mayr, Brace and Dunham postulated that there exist 
latitudinal increasing upward fluxes of ionization into protonospheric field tubes, 
which increase with latitude and which reduce the proton density thus raising the 
H + - 0 + transition level. In order to maintain these fluxes, it was concluded, that 
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ionization has to be removed out of the field tubes through diffusion across field 
lines, as the effect of a density gradient across field lines maintained through an 
escape of protons at higher latitude, where the field line configuration permits it. 
This density gradient across high latitude field lines was considered to be the 
density plateau constituting the plasmapause. Thus a consistent concept was de- 
scribed, relating plasmapause and ion compositional features at high latitudes. 

In this paper the above described concept is investigated. It will be shown 
that a certain diffusion mechanism postulated by Bohm (1949) can quantitatively 
account for the above described density plateau and its relation to the ion com- 
position at high latitudes. 

THEORY 

In the following section a theoretical formulation of our concept is presented. 
The basic equations are that of particle continuity for H + and 0 + . Because pro- 
tons populate primarily the magnetosphere their continuity equation (1) is 
discussed 


3 [ H + ] 

dt 


9 

-g K [H] 


[ 0 + ] - K [ H + ] 


[0] - div j H + - 0 


(1) 


including charge exchange between 0 + and H + as the major source of protons. 
K = is the charge exchange coefficient, Jh+ the diffusion flux of H + . 
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If diffusion across field lines is included, the diffusion term in Eq. (1) 


becomes 


div j H + - V, j x + V„ j„ (2) 

where, omitting the subscript H + , the subscripts 1 and II indicate flux com- 
ponents and derivatives perpendicular and parallel to field lines. j ,_ and j„ 
are related through motion equations, to the density gradients across, respec- 
tively parallel field lines. Thus Eq. (1) is a partial second order differential 
equation describing the two dimensional proton distribution along and across 
field lines. 

As boundary conditions we can assume that the diffusion flux and density 
gradient along the field line is zero at the equator implying an equinox condition. 
Through the charge exchange reaction, [H + ] is given by an algebraic relation to 
[0 + ] in the charge exchange region at low altitudes, thus providing the second 
boundary condition for the distribution along field lines. 

For the distribution across field lines we can assume that at a field line 
along which an escape of ionization is possible, the density distribution is 
very low. Along some inner field line we can assume that the diffusion 
flux across it is negligibly small. Thus, theoretically the problem is formulated. 
A rigorous treatment would involve very extensive numerical calculations, so 
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a simplifying approach is made with the aim to get a quantitative insight into th< 
importance of diffusion across field lines. 

Due to the lack of chemical sinks and source in the protonosphere, the con- 
tinuity equation [Eq. (1)] reduces to 


div Jh ~ V ± jj. + j„ = 0 


( 4 ) 


Without considering the divergence of field tubes Eq. (4) becomes 


_a_ . b 

dz ^ + 3 s j" 


0 




where s is the distance of field lines along and z is the distance of lines across 
the magnetic field. 
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Furthermore we assume that ~~ j x is constant along s because of the 
large scaleheightof the protons. As a consequence of this simplifying assump- 
tion Eq. (5) can easily be integrated over s, thus leading to 


(s e “ S ) 77 i-L - ill (s) 


(0) 


when s e is the distance to the equator where j„ (s e ) = 0. 

For jj_ we have employed two different formulas, each of them describing 
different kinds of diffusion mechanisms. Their implications are discussed. 
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1. Regular Diffusion 


With this notation we understand the drift of plasma across field lines due 
to collisions between electrons and ions. Referring to Spitzer (1956) we find 


1.29 x 10 13 InAN 3 


T 3/2 B 2 

e 


Bl p 


when neglecting the effect of the gravitational acceleration which seems justified 
for higher L values. Employing for p the pressure of the electron-ion gas 
P = 2kTN (N is electron concentration) and assuming InA = 16 and T = 5 x 10 3 , 
yields 


C A A on 

j . = -4 x 10 56 z 6 N -r— 

o z 


where we adopted 


B n z 3 

0 e 


with B 0 = 0.3 gauss and z the earth radius. 


Employing Eq. (8) in Eq. (6) leads to 


2.5 x 10 55 

Jll V 

(s “ s) 
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Up to this point we did not pay any attention to j„ , the ionization flux 
parallel to the magnetic field line. It is essential for our argument to discuss 
this quantity in some detail. 

The upward flux of ionization into the protonosphere, the flux j„ , can de- 
pending on its magnitude, very strongly decrease the proton concentration above 
the charge exchange region. This is the case if the amount of protons flowing 
upwards is comparable to the net production rate of protons in the charge ex- 
change region. Geisler and Bowhill (1965) have described this dependence in 
some detail for the proton density in the diffusion barrier. At this level of our 
investigation we limit our argumentation by stating that the proton density in the 
protonosphere is generally speaking 

N - N (j„, [0], [H], [0 + ], K) (10) 

indicating its dependence on neutral composition, oxygen ion density, charge ex- 
change rate coefficient and above all, indicating its dependence on the transport 
flux j„ . By concentrating our attention on the dependence between j„ and N we 
can write 


j ii 


J„ (N). 


( 11 ) 


Going with (11) into (9) yields 
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Equation (12) describes the coupling between the density variation across 
the magnetic field lines and the density on the field line. 

To put it in a more illustrative and physical picture one could interpret (12) 
by saying that the net diffusion flux out and across a field tube according to a 
density variation, has to be balanced by an upward flux of ions into the tube, 
which in turn effects the density in the tube itself. 

We have solved Eq. (12) by assuming that within the range of our solution 
j it is constant. This assumption is justified by the fact, that according to Geisler 
and Bowhill (1965), the density is a very sensitive function of j„ in the vicinity 
of the critical flux. Thus as a consequence of this assumption, our results will 
be valid only for fluxes which are in the order of the critical flux. This essen- 
tial limitation is to be stressed. The implications are of importance for a later 
argument. 

With the above assumption the integration of Eq. (12) leads to 
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Nj is assumed as the proton density at high latitudes where an escape of ioniza- 
tion is possible. 

A very simple estimation of the escape effect may help to determine an ap- 
proximate value of Nj. Following Jastrow and Rasool (1965) we find that the 
density at the escape level is in the order of 10 s /cc. Thus we learn, that the 
escape of H + becomes effective immediately above the charge exchange region 
and so we may adopt for the escape level (somewhere between 500 and 1000 km) 
a temperature of 3.10 3o K. This leads according to Jastrow and Rasool to the 
escape flux 

j e =.3X10 S • N t 

which is in principle equivalent to the previously defined flux j „ . 

A very simple argument allows to evaluate the above relation. Obviously, 
the upper limit of N x is determined by the critical flux, which represents an 
upper limit of the proton flux. So we find, that for an assumed critical flux of 
10 6 /cm 2 sec, Nj is in the order of 10/cc, which is an extremely low value. On 
the other side, if we allow the flux j e to be appreciably smaller than the critical 
flux, it would not effect the H + density and thus N x would be a sort of diffusive 
equilibrium value which is at the escape level much larger than the permitted 
maximum value of 10/cc. 
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From this we conclude that j e has to be assumed to be in the vicinity of the 
critical flux and thus N x may be as low as 10/cc due to escape of ionization. 

Returning to Eq. (13) we can now define our boundary conditions. The fol- 
lowing boundary conditions were employed which are considered to represent 
typical features of the plasmapause. 

/ BN 2 \ 

f ° r Z 0 = 4r e > ITT/ = 0 

(14) 

and for z x = 6r e , N x = 10/cc 


This means that we assume, that at L = 6 escape of ionization reduces the 
proton density to 10/cc and that below L = 4 this escape of ionization is almost 
ineffective thus reducing the diffusion flux across field lines to a negligible 
small amount. 

With Eq. (14), Eq. (13) becomes 


N 2 


100 + j„ 


2.5 x 10 55 1 

(s e - s ) 10 



(5 z - 4z 0 ) 




• (15) 


which can be interpreted as a relation between n and j l( . 

Postulating that according to satellite observations N (z Q ) is in the order of 
10 3 /cc, we get a density distribution as shown in Fig. 1 (dashed line), that has 


9 



the shape of distributions shown from OGOA satellite. However it is essential, 
that the flux j ,, has to be assumed in the order of 2/cm 2 sec, which is by several 
orders of magnitudes smaller than the critical flux, lieing in the range from 10 s 
to 10 8 /cm 2 sec -1 . Such a small flux, being only a fraction of the critical flux has 
a negligible effect on the proton density. 

The quite reasonable density decrease shown in Fig. 1 causes such a small 
diffusion flux across field lines, that the balancing proton flux into the field tubes 
is ineffective in reducing the density consistently with the density decrease 
across field lines. The diffusion mechanism employed here is inappropriate in 
describing the plasmapause. 

2. Irregular Diffusion 

Because of this failure, the question arises whether another diffusion 
mechanism could possibly be responsible for the observed density distribution. 

Spitzer (1956) also refers to a different diffusion mechanism, which was 
suggested by Bohm (1949). This diffusion is produced by a sort of turbulence or 
plasma oscillations with randomly varying electrostatic fields. Braginskii (1965) 
points out, that when cot » 1, plasma turbulence can in principle cause a strong- 
enhancement of perpendicular (to the magnetic field) transport processes. 
Though little knowledge exists on the conditions under which plasma oscillations 
are excited, we have employed Bohm's diffusion coefficient. As will be shown, it 
offers an appropriate mean in describing features of the plasmapause. 
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With Bohm's diffusion coefficient the flux perpendicular to field lines is 


3.9 x 10 18 Bp 
B TT ' 


( 16 ) 


Employing formula (16) instead of (7), we proceed like in the previous sec- 
tion, adopting again the boundary conditions (14). The result is the relation 


N 


10 + 


j„ 1.4 x 10 +19 
(s - s) 




• ( 1 ?) 


Postulating again, that at L = 4,N = 10 3 /cc, we get a distribution also shown 
in Fig. 1 (solid line). The main difference to the preceding result is, that now j M 
is in the order of 10 6 /cm sec. This is a flux which is in the vicinity of the 
critical flux, thus it can effectively decrease the proton density and thus the 
density variation shown in Fig. 1 (solid line) can be consistent to the density 
variation along the corresponding field lines, which was not the case employing 
the regular diffusion mechanism. 

As pointed out previously, Eq. (15) and Eq. (17) contain the assumption, that 
the flux j„ is constant. This is valid if j„ is close to the critical flux, because 
then the proton concentration is a very sensitive function of j M and thus, in turn, 
this flux is n early constant for a wide range of [H + ] . Obviously this condition 
was not fulfilled employing the regular diffusion mechanism with (15) and conse- 
quently the dashed distribution in Fig. 1 is invalid.With Bohm's diffusion 
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coefficient in (17) however the resulting flux j M comes into the range of the 
critical flux and thus can satisfy the above condition. Therefore the solid dis- 
tribution in Fig. 1 may be considered as reasonable. 

Comparing our calculated density variation with those seen from OGO-A 
(Taylor et al., 1965) shows relatively good agreement. However our result ap- 
pears to be in contrast to IMP 11 results by Serbu and Maier (1966) which 
indicate a density decrease with z -3 . 

In addition, the decrease of [h + ] due to escape and diffusion across field 
lines increases the 0 + - H + transition level, thus providing a possible explana- 
tion for the discrepancy between the cited ion composition measurements 
(Taylor, Thomas) and theoretical results (Mayr et al.) on the ion composition at 
higher latitudes. 

The conditions under which Bohm's diffusion mechanism is appropriate are 
not well understood. It seems plausible, that under magnetospheric conditions 
of low density and magnetic field strength, plasma instabilities and oscillations 
may relatively easily excited, so for example by the interaction with electron 
beams. Yet, the applicability of this diffusion mechanism in the magnetosphere 
is not established. Our phenomenological approach merely points at the possi- 
bility that it might be valid there. 
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Figure 1. Density distribution between L values 4 and 6 tor different diffusion mechanisms. Note the very different values of the 
proton flux jj[, on which account clearly the irregular diffusion has to be considered as appropriate. 



